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ABSTRACT 

The  axisymmetric  elastic-plastic  torsion  of  a shaft  subject  to  the  von  Mises 
yield  criterion  is  considered.  The  problem  is  reformulated  as  a variational 
inequality  and  it  is  proved  that  the  problem  has  a unique  solution.  Some  proper- 
ties of  the  solution  are  derived. 
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SIGNIFICANCE  AND  EXPLANATION 

When  an  axially  symmetric  shaft  is  subjected  to  a small  torque  the  shaft 
deforms  elastically.  As  the  torque  increases,  the  maximum  stress  reaches  the 
largest  value  permissible  in  elastic  deformation  and  a plastic  enclave  forms 
and  grows. 

In  this  paper  the  problem  of  the  elastic-plastic  torsion  of  a shaft  is  re- 
formulated as  a variational  inequality.  This  is  mathematically  equivalent  to  the 
principle  of  Haar  and  von  Karman  according  to  which  the  strain  energy  must  be 
minimized  subject  to  the  constraint  that  the  stress  should  not  exceed  its  per- 
missible limit. 

The  advantages  of  formulating  the  problem  as  a variational  inequality  are: 

(i)  The  elastic  and  plastic  regions  are  treated  in  a unified  manner  and  there 
is  no  noed  to  determine  the  boundary  of  the  plastic  region. 

(ii)  Mathematical  questions,  such  as  existence  and  uniqueness,  are  readily 
answered . 

(iii)  The  variational  inequality  lends  itself  to  numerical  approximation. 

Wo  establish  existence  and  uniqueness  of  the  solution,  and  also  obtain  bounds 
for  the  size  of  the  plastic  region.  Numerical  results  for  a two-diameter  shaft 
will  bo  given  in  a later  paper. 
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THE  SOLUTION  OF  THE  AXISYMMETRIC  ELASTIC-PLASTIC  TORSION 
OF  A SHAFT  USING  VARIATIONAL  INEQUALITIES 


* 

C.  W.  Cryer 


1.  Classical  formulation  of  the  problem. 

The  problem  to  be  considered  is  shown  in  Figure  1.1.  Equal  and  opposite  toroues  T are 
applied  to  the  ends  of  a shaft  of  length  L which  is  axially  symmetric  about  the  x^-axis 
and  has  (variable)  radius  Rtx^ . 


Figure  1.1:  A circular  shaft  of  varying  diameter. 


Because  of  axial  symmetry  it  suffices  to  consider  the  problem  in  the  two-dimensional 
domain 

0 = {x  = (x^x^  ! 0 £ < L;  0 < x2  £ Rtx^  ) , (1.1) 

corresponding  to  the  cross-section  of  the  shaft. 

The  boundary  r of  0 consists  of  three  parts:  r^,  and  u r22  as  shown 

in  Figure  1.2.  T21  and  I’22  are  parallel  to  the  x.,-axis.  is  the  curve  x2  - R(x^) , 

0 £ x^  £ L.  rQ  is  a sequent  of  the  x^-axis. 

As  regards  the  boundary  it  is  assumed  that: 

2 

(i)  R « C (0/L) , that  is,  R is  twice  continuously  differentiable.  This  assumption 
allows  us  to  prove  that  the  solution  is  differentiable  (see  Theorem  5.9). 
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(ii)  dR/dx^  = d R/dx.,  = 0 when  x^  = 0 and  x^  = L.  This  is  true  if  the  shaft  has  con- 
stant radius  near  its  ends  as  often  hc.-oens  in  practice.  This  assumption  allows  us  to 
reflect  Q in  and  T ^ and  obtain  a smooth  solution  in  the  enlarged  domain  (see 

Lemma  5.1) . 

.iii)  dR/dx^  >_  0,  so  that  is  of  the  form  shown  in  Figure  1.2.  This  assumption  allows  us 

to  conclude  that  R(x^)  >_  R(0)  for  x^  c [0,L).  It  also  allows  us  to  conclude  that 
only  one  characteristic  passes  through  each  point  on  T21  (see  Theorem  4.2). 
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Figure  1.2:  Cross-section  of  an  axisymnetric  shaft. 


In  analogy  with  the  theory  of  torsion  of  prismatic  bars  due  to  Saint-Venant  (Love  (1944, 
p.  311)),  it.  is  assumed  that  the  only  non-zero  stresses  are  shear  stresses  on  the  planes  Q. 

It  can  then  be  shown  (Love  (1944,  p.  325),  Eddy  and  Shaw  (1949),  Zienkiewicz  and  Cheung  [1967]) 
that  the  problem  reduces  to  finding  a stress  function  u.  The  stress  components 
and  x . 'x  are  given  in  terms  of  u by 


T23  = ’ u'l/(x2)  ' 
t13  = + u'2/(x2)2  ' 
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re  u,j  » 3u/3Xj,  The  stress  q is  given  by 


2,1/2.  „ 1/2 
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(x2)*  (x2) 


(1.3) 


When  the  torque  3'  is  small , the  stresses  are  small  and  the  response  is  elastic.  As  1' 
increases  a small  plastic  enclave  forms.  In  general , 0 is  divided  into  two  subregions , the 

elastic  region  dfi  and  the  plastic  region  ft  . The  unknown  free  boundary  between  1)^  ami 
1)  la  denoted  by  1'^  (see  Figure  1.2). 

In  ti  the  material  is  elastic  and  u satisfies  the  differential  equation 


An 


■ - (u,  /(X  )3),  1 - f rj-f— -T^rY  in  1) 

1 2 * 1-1  8xi  \(x2)3  3xi/  * 


(1.4) 


The  matoiial  la  assumed  to  yield  according  to  the  criterion  of  von  Mines)  that  is,  the 
material  yields  when  the  aliens  q reaches  the  maximum  permissible  value  k (a  qivon  con- 
stant). Thus, 


|qrad  u|  ^ kx.,,  in  11^  , 


2 

I grad  u|  - kx^,  in  11  . 

The  Innindaty  conditions  foi  u on  I'  ate  (see  Figure  1.2)  1 

u - 0,  on  1' 


0 


u - T/2«  on 


Jn  u„  w*i  ‘ 0,1  ri 


( 1 . r>) 

(1.0) 

(1.7) 
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Condition  ti."M  .11  isos  fiom  the  axial  symmotiy  of  the  ptoblem.  Pond  it  ion  tl.tl)  evptussr.s  the 
tact  that  the  total  torque  is  T and  that  tluue  is  no  1 1 act  ion  on  the  outei  surface  1 ^ . 
Condition  t 1 . Ui  expresses  the  assumption  that  at  the  ends  of  the  shaft  the  st  teases  cot  1 ott|suvl 
to  a puie  1 01  quo  so  that  i,(  - u, ^ - 0. 
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Tho  formulation  of  the  problom  is  completed  by  the  requirement  that  u and  its  first 
derivatives  bo  continuous  across  T^.  The  problem  defined  by  (1.4)  to  (1.9)  will  bo  called 
the  Classical  Problem. 

The  remainder  of  the  paper  is  organized  as  follows.  In  the  remainder  of  this  section  we 
make  some  brief  remarks  about  related  work  in  the  literature,  and  indicate  the  reasons  for 
choosing  the  method  of  analysis  used  in  this  paper.  In  section  2 we  introduce  certain  weighted 
Soiolov  spaces)  in  soction  3 wo  analyse  the  one-dimensional  problem;  in  section  4 the  classi- 
cal problem  i3  reformulated  as  a variational  inequality;  and  in  soctions  5 and  6 the  existence 
of  a solution  and  various  properties  thereof  are  proved. 

Numerical  results  will  appear  in  a subsequent  paper  (Cryer  11979a) ). 

In  recent  years  the  elastic-plastic  torsion  of  cylindrical  bars  has  been  intensively 
studied;  see  Ting  119731;  hanehon  11974);  for  other  references  see  Cryer  11977,  soction 
I.S.3.1).  if  the  cross-section  of  the  bar  is  denoted  by  it,  then  it  is  required  to  find  a 
stross  function  4 such  that 


A|  - - 


11 


<•>*> 


+20-0,  in  fi 


| grad  $ | » k,  in  0 


(1.10) 


$ - 0,  on  3i)  . 

Here,  the  constant  k denotes  the  maximum  stress,  and  the  constant  0 s 0 denotes  the  angle 
of  twist  per  unit  length  of  the  bar,  while  ilp  and  denote  the  plastic  and  elastic  re- 

gions, respectively. 

There  are  close  sim'larities  between  the  problem  considered  in  this  paper  and  the  pioblom 
(1.10),  but  there  are  al.o  two  important  differences; 

(i)  The  differential  o;  orator  A of  (1.4)  can  be  written  in  several  forms. 

-3 


Au  “ div(x_  grad  u) 


- X.  , Au 


2 n 
3 u ^ 9 **u 


3x' 


-L  J?iL 

3x_ 


- x,  Au  = x„ 
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(1.11) 


(1.12' 


(1.11) 


•a 

.i 


(ii)  Tho  boundary  conditions  (1.7)  through  (1.9)  area  combination  of  Dirichlot  and  Neumann 
conditions  while  tho  boundary  conditions  for  (1.10)  are  Dirichlot. 

The  singularity  of  the  operator  A is  the  most  significant  difference  between  the  pre- 
sent problem  and  problem  (1.10).  There  is  an  extensive  literature  on  degenerate  elliptic 
equations  (Visik  (1964),  Oleinik  and  Radkovic  (1973),  Pichera  [1956,  1960],  Kohn  and  Nirenborg 
(1967,  1967al  , Baouondi  and  Goulaouic  (1972)). 

Unfortunately,  much  of  tho  literature  is  not  applicable  to  tho  problem  ir.  hand.  One 
reason  for  this  is  the  following.  Tho  equation  (1.13)  is  degenerate  on  However,  the 

inner  product  of  the  coefficients  of  tho  first  order  derivatives  with  the  outward  normal  on 
r is 

(0)  v (0)  + (-3)  v (-1)  - 3 , 


which  is  positive  so  that  boundary  cond'tions  must  be  imposod  on  )'  (Pichera  (I960),  Oleinik 
and  Radkovic  (1973),  Friedman  and  S'insky  (1973)).  On  tho  other  hand,  for  the  equation 


the  inner-product  of  the  first  order  coefficients  with  the  outward  normal  on  1'^  is  equal  to 
-3  so  that  no  boundary  conditions  can  be  imposod  on  1'^.  Tins  means  that  papers  on  degenerate 
elliptic  equations  in  which  only  bounds  on  the  absolute  values  of  the  coefficients  of  tire 
equation  are  imposed  (Mur  thy  and  Stampacehia  (19681,  Trudinger  (19731),  cannot  be  of  use  in 
the  present  case. 

However,  tho  operator  A gives  rise  to  generalised  axially  symmetric  potentials  which 
have  been  extensively  studied  (see  Weinstein  (1953),  Huber  (1954,  1955),  Quinn  and  Weinaeht. 
(197o),  Quinn  (1978),  and  tho  references  below).  Various  methods  have  been  used  to  studv 
boundary  value  problems  for  generalized  axially  symmetric  potentials: 

(a)  Max  i mum  1'  t i nc  t p 1 o ■ J.imet  (1967,  1968),  Part  or  [1965,  1965al  , l,o  (1973,  1976). 
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(b)  Perturbation  of  0.  The  problem  is  considered  in 

fl  » { (x,  ,x  ) e fl  : x > e)  , 
t l i i — 

and  than  the  limit  is  taken  as  e ■»  0.  (Scheehter  [I960],  Greenspan  and  Warten  [1962]). 

(c)  Weighted  Sobolev  spaces.  The  problem  is  reformulated  as  a minimization  problem  in  the 
space  of  functions  u such  that 


. i 2 2 

1 — (u#1  + u,2)dx1dx,  < 


0 x^ 

(Leventhal  [1973,  1975],  Jakovlev  [1966],  Nocas  (1967,  chapter  6]). 

For  the  present  problem  the  natural  setting  is  a weighted  Sobolev  space,  but  we  also  use 
the  maximum  principle  and  perturbation  of  0. 

Noncoercivo  variational  inequalities  have  been  considered  by  Lions  and  Stampacchia 
[1965],  Levy  and  Stampacchia  [1971],  and  Deuel  and  Hess  [1974],  but  none  of  these  results  are 
applicable  to  the  problem  considered  hero. 
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2.  Some  weighted  Sobolov  spaces. 

3 

Because  of  the  fceim  l/(xj  in  the  operator  A defined  by  (1.4),  it  is  necessary  to 
introduce  Sobolov  spaces  with  a weight  function 

P ( x)  « pfx^.Xj)  « (x2)'3  . (2.1) 

There  is  an  extensive  literature  on  wel <hted  Sobolov  spaces  (Necas  (1967),  Kudrjacev 
(1974),  Kadlec  (1966),  Kadlec  and  Kufner  (1966,  1967),  Kufner  (1965,  1965a,  1969),  Jakovlev 
(1966)).  The  present  problem  presents  several  aspects  which,  taken  together,  are  not  covered 
in  the  literature: 

(i)  The  weight  function  p involves  the  distanco  to  the  plane  x2  " 0. 

(il)  p « x/  whoreaa  most  references  consider  the  case  p » x”,  n ^ 0. 

(ili)  The  boundary  conditions  on  3ft  are  of  the  third  kind  (Dirichlet  and  Neumann). 

Thu  results  of  this  section  hold  whenever  !1  is  of  type  N*0*'1,  that  is,  11  is  a 
bounded  domain  whose  boundary  is  LipachitE  continuous  (Necas  (1967,  p.  55)).  This  condition 
is  satisfied  as  long  as  !‘j  consists  of  a finite  number  of  LipachitE  continuous  curvos,  with- 
out cusps,  and  is  certainly  satisfied  when  0 is  as  in  Figure  1.2. 

t,"(ft)  and  wm,l>(ft)  donote  the  usual  Lebesgue  spaces  and  Sobolev  spaces  defined  over  si. 

2 

We  denote  by  L * L“(sl)  the  real  linear  space  of  real  measurable  function'  v defined 
on  s)  with  finite  norm 


! 

i 

I 

I 

» 

| 

! 

i 

4 

! 

j 

\ 

i 

* 

! 

l 

4 


1/2  1/2 

I v ; L 1 1 " ||v;I."(sl)  II  ■ (/  p(x)v*  dx)  » | / x~3v“  dx) 

P ft  ft  " 


(2.2) 


1/  ’ 2 

Thus,  v < L iff  p v , L (ft) . We  assert  that  L is  complete.  To  see  this  lot 


‘V 

bo  a Cauchy  sequence  in 

L.  Then  (p 

1/2  i , 

v 1 n 

n 

i a Cauchy  sequence 

2 

in  l/(Sl).  Since 

i 

4 

! 

! 

I '(ft) 

1/2 

is  v’omplete,  p v » 

u in  l.2(ft) 

for  some 

U l 1,'  (ft) . 

Thus 

1/2  1/2 
a v ► p v in 
n 

! 2 (ft) 

-1/2 

wlieie  v » up  < L. 

That  is,  v 

n 

► v in 

L,  so  that 

L is 

indeed  complete. 

! 

1 2 

We  denote  by  W « K * “(iJ)  the  spare  of  functions  v < L with  generalized  derivatives 
P 

fx  0 v,  1 v \ v 2/  which  also  belong  to  L.  As  noun  we  take 


I v ; W j | - ||  vsW*  ’ * (ft)  ||  » l II v ; 1. 1 1 * ■»  l !|d.v;l||  *") 
1 1=1 


, 1/2 


(2.3) 
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We  assert  that  W 13  a Banach  space.  To  see  this,  let  (v^J  be  a Cauchy  sequence  in  W. 

1/2  1/2  2 

Thon,  by  the  arguments  of  the  previous  paragraph,  p v **  P w c L (tt)  for  some  w t L, 

I/-)  1/*^  2 

while  p D.v  p w.  < L (fl)  for  some  w.  e L,  1 < i < 2.  We  must  show  that 
ini  i — — 

w.  «•  O^w.  To  do  so,  choose  a teat  function  s?  e 0( fl) , the  sot  of  infinitely  differentiable 
functions  with  compact  support  in  fl.  By  definition, 

/ (D.v  )\?dx  » - / v (D,  'f) dx  . 

0 1 n 0 n 1 

Now  v*  has  compact  support  on  fl.  Thus  *>  0 outside  some  contact  subsot  fl_  of  fl.  On 

1/2  1/2  "> 

fl  wo  have  that,  for  some  t,  x,  > < > 0.  Since  p v ■»  p w in  L“(fl),  we  conclude 

c 2 — n 

that  v |fl  • w|fl  in  L*(fl  ).  Similarly,  D.v  |fl  •»  w.  |fl  in  L‘(fl  ).  Thus, 
n c c c inc  1 c c 


/ w(D.  sf)dx  <■  / w(D.*)dx  » lim  / v (D.  tf)dx 

* rt  ‘ _ n • 


n*"*  !) 


■ • Urn  | (D.v  Vdx  • - / w.  v>dx 
x n a * 


n +»  U 


- - / w . *>  dx  , 
fl  1 

and  wo  conclude  that  indeed  w^  » D.w. 

The  preceding  arguments  used  only  the  fact  that  p is  continuous  and  positive  in  .). 
The  arguments  which  follow  use  the  fact  that  p “ x/. 

We  denote  by  V « V^'*(fl)  the  set  of  real  measurable  functions  v defined  on  fl  such 
that  x^v  « L and  v has  weak  derivatives  D.v  c l..  As  norm,  we  take 


1 v ; V 1 1 « ||v)Vl,2(fl)  ||  = [||x~lv;U|j2  +■  l ||d.  V.-Lll  2) 

i»l 


Using  the  arguments  previously  applied  to  W it  follows  that  V is  a Banach  space. 
If  v , V then  v « W and 

1| v; W 11  v_  (max  (l  *■  x,)  I ||v;v|| 

fl 

so  that  V ca;.  be  imbedded  in  W. 


! 

f 


For  small  positive  h let  be  the  strip 

■ {x  t f!  : 0 < x2  < h)  » (2.6) 

and  let 

nh  » {x  c n : x2  > h)  - fl/Sh  . (2.7) 

Let  Cq(T^)  » CQ(R”|n)rj)  denote  the  set  of  restrictions  to  0 of  functions  which  aro 

2 

infinitely  differentiable  in  R and  which  vanish  in  some  neighborhood  of  F^.  In  particular, 
if  ^ £ Cq { Tq)  then  <f  vanishes  in  some  S^.  We  denote  by  ^W  » *\(^'2(fi)  the  completion 
in  W of  C*(ro),and  sot 

II w « °W||  “ 1 l lib.viLlI2)  ' . (2.8) 

i-1 


Theorems  2.2  and  2.3  below  aro  based  on  rosults  due  to  Kadloc  and  Kufnor  (1966). 

Wo  use  the  following  inequality  duo  to  Hardy  (Hardy,  Littlowood,  and  Polya  1193-1,  p.245)). 

1-omma  2.1.  (Hardy) 

If  p>l,  n<p-l,  and  g(t)  is  a measurable  function  on  (0,«>)  such  that 


/ |g(t) |pt°dt  < ® , 

0 


then 


ft 

/ |g(s)|ds 
0 


ta'pdt 


* NT 


|g(t)|Vdt 


n 


Theorem  2.2. 


The  norms 


V = W . 

||  w ; V 1 1 ~ = ! o(x1'w“'  + Jgrad  w|“)dx  , 

si 


||  w;  W || 


0 


/ p[w‘  + [grad  w|2)<ix  , 

s) 
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|w)W°||  «>  / pjgrad  w|  dx  , 

fi 


are  equivalent  on  W,  and  satisfy 


T llw>vl|2  £ || w» °wj| 2 £ ||wjw|| 2 £ max  (1  + x )2||w)V||2  , 


If  w c W thon  wfx^x^  •*•  0 as  x2  ■*•  0 for  almost  all  x^.  Indeed, 


w(x1#x2)|  £ -j=  / ~ |d2w(x1#b)  l^ds  , a.e 

0 x_ 


2 . 1/2  5/2 


/ -r  w2dx  £ i ||  wj  °W| 
i!  x_ 


Proof ; Lot  w < W.  Thon  w belongs  to  tho  Sobolev  space  H (fl)  and  so  wlx^*)  is 
absolutely  continuous  as  a function  of  x2  for  almost  all  x^  {Morroy  [1966,  p.  66]) 


wfXj.t)  - w (x^ , s)  - / D2w(xl(u)du  . 


Furthormoro,  since  ||wiw||  < »>  it  follows  from  Fubini's  theorem  that 


1 > 2 
/ — |d2w(xi(u) | du 


for  almost  all  x^.  Thus,  using  Holder's  inequality, 


3/2  1. 1 3/2 


Iwtx^.t)  * w(x^,s)|  <f  |o2w(x1,u) I . |u|  .|ul 


rt,  ai/ir  t 

ij/  -y  |d2w(xi(u)  |“du  j / |u|  du 

l_s  u J vs 

, I X2  . , "j  1/2  . 

<-=\(  — |dow(x,,u)|  du  |t-s|  / 


««-»-»_ - — ^^Mmasea ww?»a>  1- — — 


from  which  we  conclude  chat  w(x,,0)  = lim  wfx^s)  exists  for  almost  all  x^  However,  from 

s+0 


Fubini's  theorem, 


R(x2) 


/ — |w(x^,u) |^du  < “ , a.e.  , 


so  that  w(Xj,0)  ■ 0 a.e.  Indeed,  wo  have  (2.10). 


Applying  Lemma  2.1  with  g » DjW,  a - -3,  and  p - 2,  we  see  that,  for  almost  all  xJ( 


Rtx^) 


Rfx^ 


/ \ |w(x1(x2) |2dx2  - / 


/ D?w(x.,8)ds 

0 


dx2  , 


if 


R (xx> 


'2  l 


/ jDjWtXj.s) |ds 


dx2  , 


R(x1) 


<jf  |d2w(xi(x2) |2dx2 


Integrating  with  roapoct  to  x^  wo  obtain  (2.11).  The  romaindor  of  the  Theorom  now  follows 
immediately. 

□ 


Remark  2.1. 

If  v » 0 on  DO  then  inequality  (2.11)  is  related  to  tho  Poincare  inequality.  For 
genoral  mixod  boundary  conditions,  one  obtains  an  inequality  such  as  (2.11)  only  when  Si 
satisfies  certain  restrictions  (Stampacchia  11969,  p.  145)). 

D 

a> 

Theorem  2.3.  Given  v t V and  c > 0 there  exists  t C^U’^)  such  that 

||  v - I>  iW  II  < t . 

CO 

If,  in  addition,  v - y c C^U’^)  for  some  constant  y then  ^ can  be  chosen  so  that 
1>  - Y < c'dy  . 

Consequently,  V « W°  = W. 
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Proof i 


W 1 

n I is*  \ 


Choose  f(t)  <■  C (H  ) such  time  0 <_f(t)  < 1,  f « 0 'or  t <_  1,  and  f “ 1 for 
t ^ 2.  Such  an  f can  ho  constructed  using  molllflors.  Lot  c - max  |f'|. 

I.ot  Kj^x^iXj)  » f(x2/h),  ao  that  0 < < 1, 


Ph(Xl'X2* 


1 0,  if  x2  < h , 

11,  if  x,  > 2h  . 


T)\on  Ph  t CQ{\'Q)  and 


\nL  Ph|  < oh"' 


For  any  v c v lot  v.  - K.  v.  Thom  v,  « Vi  v.  (x)  » 0 for  x e S.  i v(x)  » v.  (x) 
n n ti  n li  h 

for  x < Sl2hi  |v(x)  - v)((x)  | < |v(x)  | . Mao, 

K<v'Vl2  < il»  - vdiv)  + lv  ‘V’iJ'2  * 

< 2l|DiV|2  t V2|0iPh|2l  , 

ao  that 

|0i<v-vh)|2  - Ivl2  • ln  a»,  - 

|0jlv-vh>|2  < 21 |oAv|2  ♦ v2o2/h2l  , 

< 2l|t>1v|2  ,4vV/k*i.  in  s2)i/r»li  , 

Thua,  romomhorlmj  that  v » v.  • 0 in  Si„,  , 

»»  <*u 

II  (v-vh)  |W||2  < / | v 1 2dx  t 


b3h  *2 


1 l 2/ 

l-l  S 


2h 


~ |o,v|“  t <lc2  | v | 2 


dx 


Since  v t v each  integral  ia  convergent.  Since  tho  measure  of  S2^  ► 0 as  h > 0,  w con- 
clude that 

II  v“v|(iw||  > 0,  as  h * 0 . 
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1 ,,, , ,.....„  , _ f 


l*WHPPHlW‘l.Wl  '** 

*■  - - „ r v,  - ■ . ; _*■ , 


3.  Tho  ono-dlmcn3lonal  problem . 

It  la  instructive  to  consider  the  one-dimensional  problem  which  arises  when  tho  shaft  has 
constant  diameter.  In  thi3  case  u depends  only  upon  x2«  It  is  convenient  to  sot  x2  ••  r. 
We  normalize  u and  r so  that  tho  shaft  has  radius  1(  and  u ■ 1 on  the  outer  surface  of 


tho  shaft.  To  be  consistent  wo  should  sot  ft  » (0,L)  * (Oil)  but  wo  set  ft  * (Oil)  since 
no  confusion  can  arise.  We  look  for  a solution  for  which  ft  ■ (0(t)  and  ft  ■ (til)  for 


some  constant  r . 


Conditions  (1.4)  through  (1.9)  bocomoi 


3 i 1 3u. 

Au  " - (T  d7>  " °* 


l Du  I . 2 

laF*  ’ kr  • 


0 < r < r 


t < r < 1 , 


|3ui  _ , 2 
Irrl  < kr  , 


0 < r < 1 i 


u « Oi 


r » 0 i 


u - 1, 


r - 1 . 


Integrating  (3.1)  wo  see  that 


3u  . 3 

d7  “ 4ar  , 


0 < r < t 


for  some  constant  a.  Integrating  again  wo  obtain 


u " ar  + bi 


0 < r < r 


for  some  constant  b.  It  follows  from  (3.4)  that  b » 0 so  that 


0 < r < t . 


From  (3.2) , 


3u  . , 2 
3r  " 1 kr  ' 


t < r < 1 . 


Since  u is  required  to  be  continuously  difforontiublo  at  r « t,  tho  constants  a and 
Jtk  must  have  tho  same  signi  so  that  3u/3r  has  tho  same  sign  throughout  (0,1).  From  (3.4) 


and  (3.5)  wo  see  that  Du/3r  must  bo  positive. 
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Thus  (3.2)  becomes 


3u  , 2 
37=  kr  ' 


t < r < 1 


(3.8) 


Integrating  and  using  (3.5)  we  obtain 

u = kr3/3  + (1  - k/3) , T < r < 1 . 


(3.9) 


The  expressions  (3.7)  and  (3.9)  involve  two  unknown  constants  T and  a.  We  determine 
these  by  requiring  that  u and  u^_  be  continuous  at  r ■ t.  From  (3.6)  and  (3.8)  we  have 

(t  - 0)  » 4ai3  = (t  + 0)  ■ kT2  , 

so  that 

a * k/4x  . (3.10) 


From  (3.7)  and  (3.9)  we  have 

u(t  - 0)  = ax4  = u(x  + 0)  = kx3/3  + (1  - k/3) 

Substituting  from  (3.10)  and  re-arranging,  we  obtain 

T3  = 12 (k/3  - l)/k  . (3.11) 

The  solution  x of  (3.11)  depends  upon  the  value  of  k.  There  are  three  possibilities: 

1.  k < 3.  Then  x < 0.  Physically  this  means  that  the  torque  T is  too  great  and  no 
solution  exists. 

2 . k > 4.  Then 

x “ [4  — 12/k]1/3  > 1 . 

Physically  this  means  that  there  is  no  plastic  region,  and  the  analysis  must  be  modified. 

4 

Setting  a = 1 in  (3.7) , we  obtain  a solution  u = r of  the  elastic  problem  which 
satisfies  '.he  constraint  (3.3)  namely  ||^|  <_  kr2. 

3 . 3 < k < 4 . Then 

x = (4  - 12/k]i/3  e (0,1)  (3.12) 

and  there  is  both  an  elastic  region  = (0,x)  as  well  as  a plastic  region  (x,l). 

From  (3.7) , (3.9) , and  (3.10) , 
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i 

i 


rt»wt  »■  » » Ki  wHiOh  «Wi 


[kr4/4r. 


In  (1  , 

e 


(3.13) 


Lkr  /3  + (1  - k/3) , in  ft  . 

r 


We  now  show  that  u,  as  given  by  (3.13),  satisfies  two  alternative  formulations  of  the 
problem. 

Direct  computation  shows  that 


Au 


, in  flo  , 


lk/r  > 0,  in  . 


(3.14) 


I.et  <J)  bo  such  that 


grad  \|>  - ■ kr2, 

td)  “ 1 , 


0 < r '•  1 


(3.15) 


so  that 


<!>  - kr  /3  + (1  - k/3) , 0 < r < 1 


(3.16) 


\Ji  is  called  the  obstacle. 

3 

Noting  from  (3.11)  that  k/3  - 1 « kt  /12,  direct  computation  shows  that 
f k(r-t) 2(3r2+2rt+T2)/12T  > 0,  in  ft  - (0,r)  , 


U-\|l 


(3.17) 


0 , 


in  ft  " (t,1) 

D 


Combining  (3.14)  and  (3.17)  it  follows  that  u satisfies  the  one-dimensional  Complementary 
Problem i 


Au  >_  0,  in  ft  , 

u - <l>  j>  0,  in  ft  , 

(Au)  (u  - \)j)  = 0,  in  ft  . 


(3.18) 


Now,  with  the  notation  of  section  2 lot 


Set 


V » V1(2(fi)  <■  °W1,2(ft)  . 

P P 


K “ (v  c V s v(l)  « 1 ; v(r)  > ^(r)  for  r < (0,1)) 


(3.19) 


(1.20) 


-16- 


V c w ' (Oil)  ami  so  If  v « V than  v la  equivalent  to  an  absolutely  continuous 
function.  Thus,  statements  such  as  v(l)  « 1 in  the  definition  of  K can  be  interpreted  in 
the  classical  sense.  Furthermore,  since 

i 3. 

J —r  v dr  < **  , 

0 r 

we  see  that  the  condition 


v(0)  - 0 


is  satisfied  by  all  v e V. 


I,ot  a be  the  bilinear  function  on  V x V, 


(3.21) 


i , 

a(u,v)  « / — r u (r)v  (r)dr 
n r r 


(3.22) 


Then,  for  any  v c K,  and  remembering  that  u » ar4  in  (0,\), 

1 , 

at'WV  - u)  « / u (v  - u ) dr  , 

0 rJ  * r 1 

“ i J>  Vvv  - Vdr  f f'  ’T  VVr  - Vdr  • 


Integrating  by  parts, 


a(u,v  - u)  - 


1 lr  1 

(v  - u)  ~ u + / (v  - u) Au  dr  t 

L r ‘Jo  0 


♦ (v  - u) 


, 'll  1 

^ ur  * / 
r i 


(v  - u)Au  dr 


Since  v(0)  - u(0)  - 0,  v(l)  » u(l)  - l,  An  - 0 in  l)o  - (0,  t ) , u - in  s!^  - (t.l),  and 
is  continuous  at  r - t,  vra  obtain 


n(u,v  - u ) - / (v  - i)»)Au  dr  . 


but,  v < K so  that  v *1',  and,  by  (3.14),  au  x 0 in  rt()  - (x,l),  so  that  u is  a 


solution  of  the  ono-dimons tonal  Variai 


Inequality i Find  u t K such  that 


a(u,v  - u)  s 0,  for  all  v , K. 


O..M) 


i 


a*^!^T?»:»SK-E»<asM», 


4 . Tho  Two-dimensional  Variational,  Inequality. 

In  tho  previous  section  it;  was  shown  by  direct  computation  that  tho  solution  u of  tho 
classical  one-dimon3ionnl  olnatic-pl iStic  problem  satisfies  tho  ono-dimensional  complementarity 
problem  (3.10)  and  tho  ono-dimensional  variational  inequality  (3.23).  This  suggests  that  wo 
consider  tho  corresponding  uwo-dimansional  problems. 

Tits  two-dimensional  c wpiomontarity  problem  is  vory  useful  conceptually,  and  also  very 
helpful  when  one  consider , numerical  approximations.  However,  this  problem  gives  rise  to 
technical  difficulties  s.nco  it  is  necessary  to  carefully  define  tno  moaning  of  statements 
such  as  (Au)  (u  - <J0  >0.  This  can  be  done,  but  wo  will  not  do  so  here. 

In  contrast,  the  two-dimensional  variational  inequality  is  relatively  easy  to  apply  since 
wo  can  use  tho  following  fundamental  result  of  Stampacchia  [1964) : 

Theorom  4.1i  Lot  V be  a real  llilbort  space.  Lot  a bo  a roal  bilinoar  operator  on  V x v 
such  that  a is  coorcivo  and  continuous!  that  is,  there  are  roal  strictly  positivo  constants 
and  a2,  such  that 

■ 2 


a(v,v)  a ^ ||  v I 


for  v r,  V , 
|n(v,w)|  _<_  «_  l| v ||  || w || , for  v,w  c V . 


Lot  f bo  a roal  continuous  linear  functional  on  V.  Lot  K bo  a closed  convex  non-empty 
subsot  of  V.  Then  tho  variational  inequality:  Find  u c K such  that 

a(u,v  - u)  ^ (f,v  - u) , for  all  v c K , (4.1) 

has  a unique  solution. 

D 

General  references  on  variational  inequalities  includo:  Duvaut  and  Lions  [19721, 
Glowinski,  Lions,  and  Tromolioros  (1976),  Baiocchi  [1978],  Glowinski  [1978),  Kindorlehrer 
[1978],  and  Cryer  (1977  section  II.  11.  1979). 

In  order  to  apply  Theorom  4.1  to  tho  problem  in  hand  we  must  define  V,  a,  K and  f. 

In  doing  so,  wo  have  boon  guided  by  tho  work  of  Eddy  and  Shaw  [1949),  Brezis  and  Sibony  [1971), 
and  Lovonthal  [1973,  1975). 
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The  space  V is  taken  to  be  the  space 


V - W1,2(ft)  = °W1,2(ft)  * °W  - W 
P P 


defined  and  discussed  in  section  2.  It  was  shown  in  Theorem  2.2  that  there  are  several  equi- 


valent norms  on  V.  Here  we  use  the  norm 


l|v||  - ||v,°w)|  - [J  P | grad  v|2dx] 
ft 


The  bilinear  operator  a is  defined  on  V x V by 


a(v,w)  « / p [v i x w,x  + v,2  w,2]dx 


/ p grad  v . grad  w dx  . 
ft 


|a(v,v)  | = I|v||2  , 


a is  coercive,  and  since 


|a(v,w) I < v|  . w . 


a is  continuous. 


The  obstacle  <l>  is  the  solution  of  the  initial  value  problem  for  a first  order  partial 


differential  equation: 


| grad  t[’P  = k x2>  in  ft  , 


iji  = T/2n,  on  rxi  <_  T/2ti  in  ft  , 


where  the  restraint  <Ji  <_  T/2 it  resolves  the  ambiguity  in  the  sign  of  grad  >J>. 


The  set  K is  defined  by: 


K = {v  c v : v = T/2n  on  (in  the  sense  of  H1 (ft) ) , 


v >.  ij>  a.e.  (almost  everywhere)  in  ft) 


Here,  the  statement  'v  = T/2n  cn  in  the  sense  of  (ft)  * means  that  there  exists  a 


sequence  of  smooth  functions  {*  } such  that:  (i)  « e V;  (ii)  * - T/2n  in  a neighborhood 


of  i'x;  and  (iii) 
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|k  - v||  •*  0 as  k •►  * . 


K is  closed  and  convex. 


The  boundary  conditions  (1.7)  and  (1.9)  arc  incorporated  into  the  definition  of  Ks 
every  v e V satisfies  v *>  0 on  1'0  in  a weak  sense)  and  the  condition  3u/3n  « 0 on 
is  a ’natural’  boundary  condition  in  a variational  formulation  of  the  problem. 

Finally i the  functional  f is  aero  in  tlxs  present  problem. 

Wo  claim  that  the  Variational  Inequality  corresponding  to  the  Classical  Problem  (1.4)  - 
(1.9)  isi  Find  u c K such  that 

a(u,v  - u)  0,  for  all  v e K , (4.7) 

whore  a and  K are  as  defined  in  (4.4)  and  (4.6). 

Before  procoding  further  wo  need  some  information  about  the  function  <>. 

Theorem  4.2. 


For  x e fi, 


iMx^.Xj)  “ (kXj/3  + T/2s  - kR(0)3/3) 


For  x t t'0> 


\f> (x^ ,0)  » ll  » g(0) 


Proof i v*1  is  defined  by  (4.5).  On  1'^, 


jgrad  \J»  ] » kx“  ■ 0 , 


so  that  <>•*15  on  1'^  for  soma  constant  P. 

To  determine  6 we  note  that  <>  satisfies  the  first  order  equation 


2 2 2 4 

F(x1(s2.<>.p»q)  sp  + q - kx2«o  , 


(4.10) 


where  p « ^ and  q - <> , ^ . The  corresponding  characteristic  system  of  differential  equa- 

tions along  a trajectory  paramotorised  by  s is  (Courant  and  Hilbert  11962,  p.  }s) ) , 
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eg* 


dx1 

IT  " FP  " 2p 


dx., 

-r—  » F ■ 2q 
ds  q 1 


2Lm 


d3  - pFp  + qPq  ' 2(p2  + ^ " 2*24  1 


^ - (pF.  + P ) « 0 , 

da  * tji 


5}  - - «*F* + rx2>  - <*24 


Wo  integrate  this  system  starting  at  the  point  (0,K(0) ) where 

XjiO)  “ 0,  x,(0)  » R(0)  , p(0)  * 0,  q(0)  » kxj,  <>(d)  “ T/2#  . 


<4. .11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


From  (4.14)  we  see  that  p(s)  2 0.  It  then  follows  from  (4.11)  that  x^s)  5 0,  and  from 
(4.10)  that  q » + kx“.  Wo  aro  thus  integrating  along  1',^  and  wo  obtain  the  same  value  for 
<>  as  for  the  corresponding  one-dimensional  problem. 

By  the  appropriate  modification  of  (3.161,  we  obtain 

^(Q.xJ  » g(x0)  , (4.17) 

whore 

g(x,)  1 (kx3/3  t T/2n  - kR(0)3/3)  . (4.18) 

In  particular, 

e - ^(0.0)  » T/2n  - kR(0)J/3  . (4.18) 


It  should  bo  pointed  out  that  there  is  a hidden  complication  in  the  above  argument, 
because  if  vo  follow  the  same  approach  starting  from  the  point  (l,,  K ( l.) ) we  apparently  obtain 

$(l„0)  “ ‘172#  - KR(U3/3  t S . 

The  explanation  for  this  apparent  paradox  is  that  two  or  more  characteristics  may  intersect. 

A more  detailed  study  of  C*  (Cryer  (1979a))  shows  that  when  x,  is  small,  two  characteris- 
tics pass  through  points  (1. , x.,)  * !' 0 . This  does  not  happen  on  1\  ^ because,  as  as  readily 
seen  from  (4.11)  and  (4.14),  if,  as  in  Figure  1.2,  dR/dx0  > 0*  the  characteristics  always 
have  dXj % ^ 0 rnd  only  the  characteristic  starting  at  (0,R(0))  passes  through  the  point 
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(Q.Xj)  < ral. 


Sinco 


.'^.JL.,  — ■ ■ ■ 1. 1-,. .....I. — ■•  -■■-•■>■  “■**■*-  -5 


vc  soo  that 


Thus, 


N’»2I  £ I grad  $|  - kx2 


l^tXj.Xj)  - $(*^0)  1“  / |\|>,2|dx2  £ kx2/3 


iMx^x.,)  £ B + kXj/3  " g(x2)  • 


Q 


Remark  4.1. 


At  first  sight  it  may  socm  surprising  that  ^ is  constant  along  1'^  on  which  no  condi- 
tions wore  imposed.  This  can  bo  understood  moro  clearly  after  considering  the  dotailed  cal- 
culation of  $ as  done  by  Cryor  [1979a!. 

2 
2’ 


Alternatively,  sinco  |grad  ty|  £ kx2,  wo  know  that 


||»,°W||2  ■ / p|gra<5  $|2dx  < ~ . 

n 

It  is  known  (Kadlec  and  Kufner  [1966,  p.  469),  Lovonthal  [1973,  Lemma  6.2))  that  this  imp. '.os 
that  $ is  constant  on 

Theorem  4.3. 

Let 

kQ  - (3T/2n)R(0)’3  . (4.20) 

If  k < kQ  then  K is  empty  and  tho  variational  inequality  (4.7)  has  no  solution. 

If  k £ kg  then  (4.7)  has  a unique  solution  u. 

Proof;  If  k < kQ  then,  from  (4.8)  and  (4.9),  $ - 6 » 0 on  1'^.  Thus,  if  v e K, 

||  v ; W 1 1 2 ■>  / pv“dx  ' / p^'dx  » +">  . 

Si  SI 


Hence,  K is  empty,  and  (4.7)  has  no  solution. 
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If  k >_  kQ  then  8 < 0,  and  v * max{0,\>}  e K.  Since  K is  not  empty,  it  follows  from 
Theorem  4.1  that  the  variational  inequality  (4.7)  has  a unique  solution. 

D 

Remark  4.2. 

In  section  3 for  the  case  T « 2n  and  Rtx^)  51,  we  saw  that  there  were  three  possi- 
bilities! k < 3 (no  solution)}  3 £ k £ 4 (an  elastic-plastic  solution)}  k > 4 (an 
elastic  solution).  In  Theorem  4.3  we  only  distinguish  between  two  possibilities;  k < kQ 
(no  solution)}  k £ kQ  (either  an  elastic-plastic  solution  or  an  elastic  solution). 

Further  properties  of  the  solution  of  (4.7)  are  discussed  in  the  next  two  sections) 
these  properties  are  such  that  they  justify  our  claim  that  the  variational  inequality  (4.7)  is 
an  appropriate  extension  of  the  Classical  Problem. 


i 

I 
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5.  Regularity  of  the  solution  u of  the  variational  inequality. 
We  assume  henceforth  that 

k > ko  " (S)R<0>"3  ' 


and  set 


hQ  » lR(0> 3 - 3T/2irk  1 1//3  > 0 


(5.1) 


(5.2) 


We  prove  that  u,  tho  solution  of  the  variational  inequality  (4.7),  is  regular  by  first 
proving  that  u is  regular  in  tho  strip  near  r^,  and  then  proving  that  u is  regular 

in  ft.  , where  S.  and  ft.  are  as  in  (2.6)  and  (2.7). 

”o  h0  h0 

1 12 

We  recall  certain  properties  of  the  Sobolev  space  H (ft)  * w ' (ft)  which  are  proved, 
for  example,  by  Gilbarg  and  Trudinger  (1977,  Chapter  7 and  p.  168). 

If  v,w  e H* (ft)  then  max(v,w)  e H*(ft)  where  max(v,w)  is  defined  by 


max(v,w) (x)  ■ max(v(x) ,w(x) ) . 

If  v € H*(ft)  then,  by  definition, 

sup  v ■ inf{s  « r1  i v(x)  £ s a.e.  in  ft)  , 
ft 

sup  v = infU  t R1  i v(x)  - f £ 0 on  3ft)  . 
3ft 


where 


i(x)  - l £ 0 on  3ft  iff  max{u  - i,  0)  c 


hJ(«) 


(5.3) 


(5.4) 

(5.5) 


(5.6) 


As  a preliminary  step  in  the  analysis  we  show  that  it  is  possible  to  enlarge  the  domain 
ft  by  reflection  in  the  vertical  sides  so  as  to  avoid  the  difficulties  associated  with  r^. 
This  is  a well-known  trick  for  handling  Neumann  boundary  conditions  (see,  for  example, 
Baiocchi,  Comincioli,  Magenes  and  Pozzi  (1973),  p.  25  footnote  331).  The  arguments  are  ele- 
mentary and  rather  tedious  but  we  are  not  aware  of  any  detailed  treatment  in  the  literature. 

Let  ft0  = ft.  Let  ft*  be  the  reflection  of  ft0  in  T21  and  set  ft  = ft0  u T21  u ft* 
with  boundary  fQ  u u ?2^  u ?22  (see  Figure  5.1).  Let  u be  defined  on  ft  by  reflect- 
ion: 
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u.2<x  ,x  > » + u.  (|x1|.x  ) 


in  « . 


Tho  values  of  u,^  need  not  be  defined  on  since  it  is  a set  of  measure  zero. 


Wo  introduce  the  strips  parallel  to 


Tj  - {x  e S i jxj  <_  2d}  , 


T . *■  T,  n (3  . 

a a 


For  any  d e (0,L/2)  lot  be  a cut-off  function  with  tho  following  proportios: 

(a)  <Jd  € c“(Rn)  * 

(b)  l^d.l1  - 2d'1;  gd,2  ‘ 0 • 

(c)  gd  is  symmetric  about  x^  *>  0 , 

(d)  gd  » 1 if  |xj|  < d and  gd  » 0 if  x i Tj  . 


For  any  t C0(n) , 


1 0 

* “ *d  + *d  + *d 


0 « 0 1 OH  1 

v»d  * 9dv>»  V>d  e CQ(n  ),  and  v>d  « CQ(fl  ) , 


For  i *■  1,2,  and  u,^  defined  by  (*)  and  (**) 


/„  uti  <^d  dx  » /_  u,^  ^d  dx  0 a3  d 0 , 


bocauso  ^d  is  bounded,  u,^  <.  L (fj)  , and  tne  measure  of  i'd  goes  to  zero  as  d ■*  0.  For 
the  snmo  reason, 

i 11  V’j  , dx  “ / 11  V>,,  dx  0 as  d •*  0 . 

n *2  0 d 2 


Finally, 


4"  ^.i dx  * LJ  dx  + L « 9a#l  * dx  , 


,<1>  . T (2) 

Id  + Id  , say. 


M \ - 

As  before,  •*■  0 as  d •*■  0.  Using  the  symmetry  of  g^  and  u, 


I^2*  ” / u(x)  1<X*  f^xi,x2)  “ *t-Vx2>ldx  ^0  as  d •*•  0 

since  |gd  J < 2d"1,  and  |y)(x,,x2)  - v>(-x1»x2)J  <2d  sup  l^txll,  on  Tj. 

Thus,  for  i « 1,2,  and  utl  dofinod  by  (*)  and  (**), 

X.  - / u,.  <f  dx  » / u,,  v?  dx  + / , u,.  dx  + o(l)  , 

1 ft  1 n°  1 d n1  1 d 

“ /nu‘  i dx  + l"1**  / , u,.(-x.,x,)  V>i(x)dx  + °(1)  , 

n n 

-»  / u(i(x)  [v^C+Xj^.Xj)  + (-l)1  v^Ht^Xj)  ]dx  + 0(1)  . 

But  u,j  is  the  weak  dorivativo  of  u on  ft,  so  that 

I im  - ! “Wiv’d, l<+Xi'X2>  + V’d,i(~xi'x2)  *dx  + °(1)  ' 

" “ /<>“  *d.i  dx  - /lS  *d,i  dx  + °(1)  ' 

- - / u v>, , dx  + o(l) 
ft 

Wo  conclude  ‘•hat  the  functions  u,^  as  dofinod  by  (*)  and  (**)  aro  indeed  the  weak  derivatives 
of  ii. 

Clearly,  u c V and  |jutv||  » 21//2||u»v|| . 

Finally,  wo  note  that  if  v < K then  v^,v*  c K whore 


0 -i 

v *>  v 


ft  ' 


V1(xi,x2)  ■ V(-X1<X2)  , X € ft 


For  any  v c K, 


a(u,v-u)  • J p grad  u . grad  (v-u)dx  4 /.  p grad  u grad  (v-u)dx 
ft0  ft 
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Making  the  substitution  x^  «=  -x^  in  the  second  integral,  we  obtain, 

a(u,v-u)  ■ J p grad  u . grad(v  -u)dx  + J p grad  u . grad  (v  -u)dx  , 
ft0  ft 

- a(u,v°-u)  + atUfV^-u)  , 

> 0 > 

since  u solvos  the  variational  inequality  (4.7).  That  is,  u solves  the  variational  in- 
equality (5.0).  From  Theorom  4.1,  we  see  that  the  solution  of  tho  variational  inequality 
(5.8)  is  unique,  and  tho  lomma  follows. 

Remark  5.1. 

We  can  also  rofloct  ft  in  P22  and  obtain  results  analogous  to  thoao  of  Lemma  5.1. 
Theorom  5.2. 

u is  non-nogativo  a. a.  in  ft.  That  is, 

u - max(u,0) 

Proof:  Tho  proof  is  a modification  of  tho  proof  of  tho  weak  maximum  principle  in  Gilbarg 
and  Trudingor  [1977,  p.  168], 

Lot 

ft+  » (x  € ft  : u(x)  >_  0)  , 
ft_  *>  (x  ( ft  : u(x)  < 0)  , 
v » max(u,0) 

Then  v < tl^(ft).  Furthermore,  v « 0 and  lienee  |grad  v|  » 0 on  ft_.  Thus, 

/ p(x~2v2  + | grad  v|2]dx  - / plx^“u2  »•  |grad  u|2]dx  , 
ft  ft+ 

and  we  conclude  that  v c V. 

Similarly, 

a(v,v)  " / p|grad  vj2dx  _<  / p|grad  u|2dx  ■»  a(u,u)  . (*) 

ft  " ft+ 

Obviously,  v(x)  >_  u(x)  >_  i|i(x)  a.o. 

Wo  now  show  that  v » T/2n  on  1'^  in  the  sense  of  ll!  (ft)  . Since  u « T/2a  on  i ^ in 
tho  sense  of  M^(ft),  there  is  a sequence  (W  } with  V>  ->  u in  V,  and  'f  - T/2it  < t\U'.) 


vk  * nwtx{v>k,0)  . 


p 


lX?t 

Than  belongs  to  V and  v^  « t/2h  in  some  neighborhood  of  l'j  (bvit  v^  - T/2n  / c'”u’j)). 

Since  V>k  •>  u in  V wo  know  that  v>k(x)  * u(x)  a. a.  and  that  tho  norms  ||v>k||  arc  bounded. 

Consequently , v^ix)  v(x)  a.o.  and  tho  norms  ||v^||  nre  bounded.  V is  a Hilbert  space. 

In  a Hilbert  apace  bounded  acts  aso  woakly  sequentially  compact,  so  there  exists  a subsequence 

{v£},  which  converges  woakly  to  some  v*  e V.  Weak  convergence  in  V implies  weak  convor- 
2 

gonco  in  I.  (f))  which  in  turn  implios  pointwlse  convorgonco  a.o.  Thus,  VjMx)  -v  v'(x)  a.e., 
from  which  it  follows  that  v(x)  » v’  (x)  a.o.  and  hence  that  v <*  v' . Taking  finite  convex 
linear  combinations  of  tho  wo  obtain  a sequence  (v^)  which  converges  in  norm  to  v*  « v. 
Each  v^  is  a finite  linear  combination  of  tho  v^,  so  v^  »•  T/2ti  in  some  neighborhood  of 
rr  Finally,  applying  Tlmorom  2.3,  we  approximate  v^  by  wl\oro  i|v  e V and 

~ T/2n  f.  CQ(ri) . Since  v,  wo  conclude  tl\at  v » T/2tt  in  tho  sense  of  II1  (si). 

In  summary,  v < K. 

Now,  u solves  the  variational  inequality  (4.2),  and  so 

n(u,v-u)  « - / p | grad  u|2dx  < 0 , 

n_ 

which  implies  that:  a(u,v-u)  - 0.  Uut  then,  by  (*) , 

a(v-u.  v-u)  - e(v,v)  - a(u,u)  ••  2a(u,  v-u) 

- n(v,v)  - a(u,u)  , 

£ 0 , 

so  that  v « u. 

n 

Remark  5.2. 

barter  (1065,  p.281)  gives  an  oxampio  Involving  generalised  axially  symmetric  potentials 
whore  tho  maximum  principle  does  not  apply,  in  barter’s  oxampio,  however,  tho  region  s)  is 
symmetric  about  r and  so  tho  line  of  degeneracy  is  contained  in  d.  In  the  present  paper 
tho  line  of  degeneracy  is  on  tho  boundary  of  tt. 
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In  the  strip  S , u satisfies  the  differential  equation  Au  = 0 in  the  weak  sense; 
0 

that  is, 


Proof ; If  k > kQ  then,  from  Theorem  4.2, 


/ p grad  u . grad  v>  dx  » 0 , 

n 


(5.9) 


tytx^Xj)  <_  g(x2)  < 0 


By  Theorem  5.2  wo  know  that  u >_  0 a.o.  in  ft.  Thus,  u 0 > a.e.  in  S.  . More 

"o 

specifically,  given  a compact  subset  G of  S.  thoro  exists  c > 0 Buch  that  u >_  \|i  + c 


a.o.  in  G.  For  any  v>  c C^CG)  chooso  5 > 0 so  that  |6v>|  < t.  Thon,  v+  ™ u + Gv>  e K 
and  v » u - W t K.  Hence, 


a(u,  v -u)  - a(u,  v -u)  » f P grad  u . qrad(6v>)dx  « 0 

ft 


n 


Theorems  5.2  and  5.3  depend  on  Theorem  4.2  which  assumes  tho  specific  geometry  of 
Figure  1.2  to  evaluate  «J*.  If  1’^  is  not  as  shown  in  Figuro  1.2,  lot 

R “ lKx^)  « minllKx^)  ! 0 < < I.) 

Wo  boliovo  that  Theorem  5.3  remains  truo  if  in  the  definition  of  hQ,  U(0)  is  replaced  by 
Tho  proof  would  require  a detailed  study  of  tho  function  d>  in  tho  case  of  a general  domain, 
along  the  linos  of  tho  study  by  Ting  (I960)  for  tho  case  of  tho  torsion  of  a prismatic  bar. 
Remark  5.4. 

Theorem  5.3  provides  a bound  for  the  size  of  the  plastic  region.  This  is  particularly 
interesting  because  in  the  numerical  computations  of  F.ddy  and  Shaw  11949]  tho  plastic  region 
dips  down  near  the  corner  B on  I'  (see  Figuro  1.2),  and  it  is  far  from  clear  that  tho 
plastic  region  will  not  grow  very  rapidly  as  tho  torque  increases.  For  tho  second  problem 
considered  by  Eddy  and  Shaw  11949],  k ■»  49,  T * 6349  x 2n,  R(0)  “ 8.  In  their  numerical  cal- 
culations 
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x2  = min{x2  : (x^Xj)  £ ftp)  ^ 6*95  . 

From  (5.1),  (5.2)  and  Theorem  5.3,  kQ  ^ 37*20,  and 

x2  — h0  = 5*13  ' 

Remark  5.5. 

The  fact  that  there  is  an  elastic  strip  near  rQ  as  long  as  k > kQ  is  analogous  to 
the  situation  for  the  elastic-plastic  torsion  of  prismatic  bars  where  an  elastic  core  also 
remains  until  the  entire  bar  becomes  plastic  (Lanchon  (1974)). 

Remark  5.6. 

Our  analysis  is  not  adequate  to  handle  the  limiting  case  k * kQ.  We  conjecture  that  if 
k = kQ  then  u » iji  for  0 <.  x^  <_  x^,  where 

x^  ■ maxfx^  s R(x)  15  R(0)  for  0 < x < x^)  . 

Theorem  5.3  asserts  that  Au*=  0 in  the  weak  sense  in  S.  . We  may  thus  expect  that  u 

no 

is  regular  in  S . This  does  not  appear  to  follow  from  known  results  about  ellintic  equations, 
n0 

and  we  therefore  prove  this  by  modifying  the  corresponding  proof  for  uniformly  elliptic  equa- 
tions. We  use  Gilbarg  and  Trudinger  [1977]  as  a basic  reference,  since  this  i-s  a comprehensive 
and  readily  accessible  test. 

The  basic  idea  is  to  obtain  bounds  for  the  differences  of  the  solution  u and  then  pro- 
ceed to  the  limit. 

The  difference  quotient  in  the  x^  direction  is  defined  by: 


AjV(x)  = & v(x)  = Av(x)  = 


v(x+he,)  - v(x) 


-,  h / 0 , 


(5.10) 


where  e2  is  the  unit  vector  in  the  x^  direction.  If  v e V,  then  the  difference  quo- 
tient Ahv  is  defined  on 


(^  = f!h  = (x  £ fl  : x + te^  € 5)  for  t £ (0,hl ) 


(5.11) 


As  is  customary,  the  weak  derivatives  of  u are  denoted  by  D^u,  D^u,  etc. 
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Lemma  5.4. 


Let  v e V and  h > 0.  Let  (l1  c fij  n ft"".  Then  Av  = A^v  c I.2(ft’)  and 
1|aJv  .•  i,2(n') ||  < 11  DlV  , Lp(n)  ||  . 

Proof;  The  proof  is  a modification  of  tho  proof  of  Lemma  7.23  of  Gilbarg  and  Trudingor  (1977) 

Wo  begin  by  assuming  that  v e C^d^) . If  x e Q^,  then 

1/2  lv(*+ho  ) - v(x)) 
p1/i!(x)Av<x)  - p1/i!(x) 


h 


" r f p1/,2(x)D1v(x+to.)dt  , 

h Q 1 


so  that,  using  tho  Cauchy-Schwarz  inequality, 


p(x)  (Av(x)  I*2  < r / p(x) (D.v(x+to.)]“dt  . 

- n o 


Since  ptx+to^)  » p(x)  - x2  , 


|Av  i i.2(ft’)  ||2 


/ p(x) (Av(x) )“dx 
ft' 


< r f dx(J  p (x+to  ) (D  v(x+to. ) ] “dtl 
" ft'  0 

i h 2 

» - / dt  / pfx+te  ) (av(x+to  ))  dx 
h0ft' 

h 


r / ||d  v , L2(ft)||2dt  , 

o 


and  tho  lemma  follows  for  v e C (1^) 

But,  by  Theorem  2.3  there  exists  ^ t c“(I‘0)  such  that  ||v-<i>  j w||  <_  t for  any  c ■»  0. 

Thus , 
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|Av  , l»ptft*)  || 

< ||  w ) Lptn')||  ♦-  2 «/h  , 


< ||d,*  j ii^tn>||  -*•  2e/i)  . 

*“  X p 

< |)dxv  I Lp(rt)||  + 2</h  + < . 


Letting  < •*■  0 the  lemma  follows, 


a 


Using  the  arguments  used  to  prove  Lemma  7.24  of  Gilbarg  and  Trudinger  [1977)  we  obtain 
Lemma  5,5. 

2 

Let  v c I.  (0)  and  let  O'  cc  ft.  If  for  h < distance  (ft't  3ft)  we  have 


lldjv  j L*(ft')||  < c , 


i**mx*x 


then  v has  a weak  derivative  DjV  which  satisfies 


|Dxv  , L‘(ft)|l  < o . 


a 


Theorem  5.6. 

Lot  u c V satisfy  Au  - 0 in  S in  the  weak  sense.  Lot 

l0 

S « {x  « (x^.Xj)  ; I./4  < x^  < 3Ia/4  and  0 < x2  < hQ/2)  . 

Thom 

(i)  u|s  e (MS)  . 


(ii)  u.jjs  and  u>12|s  belong  to  L^(S), 


(Ui)  l t U'22  <1X 

5 Z 


(iv)  n < 11  (S) . u can  be  extended  as  a continuous  function  to  S,  and  u ■ 0 on 
3S  n rQ. 

4 — 

(v)  u « x2  v,  where  v is  analytic  in  S. 
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Proof i Tho  proof  is*  a modification  of  the  proofs  of  Theorem  8.0  and  8.9  of  Gllbarg  and 
Trudiisgor  [19771 . 


Proof  of  (i)  i Statement  (1)  follows  from  Corollary  8.11  of  Gilbarg  ami  Trudinger  119771 


Proof  of  (11)  * Wo  denote  by  C (3S.  /!')  tho  sot  of  functions!  which  aro  continuously  differ- 
'd 


ontiablo  In  si,  vanish  outside  s , and  vanish  in  some  neighborhood  of  as.  /!'. 

'd  0 


hot  n t C (58  /r  ) bo  such  that 

"o  0 


(a)  n <•  1 for  x « Si 

(b)  lojisl , |bj,n|  ^ Cj  for  some  constant  c^.  n in  readily  constructed  as  tho  product  of 
two  ono-dimonnional  out-off  functions. 

Par  small  positive  h.  siot 

v » n’A^u  * g*Au  , 


N“  “ I IlnD.Au  i i.^lsl)  ||' 
i-1  1 1 


Than,  for  h < dlat(»upp  n.  5S)(  /l'Q) 


N*  - J / o(b.Au)  (n'D.Auldx 
i-1  0 


“If  AlnDjU)  (D  - 1 OjUiVulds 


i-l  >1 


“ f f l (nP,u)i\”*vo  -v  * ipnbjUiAD^u) Au)dx 
i-l  si 


(*) 


where  we  have  nsed  the  identity,  valid  for  any  f t 0* (?S.  /V  ) , g t V,  and  sufficiently 

u0  U 


assail  is, 


/ (Ajf)«j  dx  - - / fA~*'g  six  . 


Since  Au  - 0 weakly  in  si.  , we  know  tls:»t  for  assy  V « 0„  (S,  ) , 

\s  1 \l 


} f t'  (bjU)b.  >>  six  - 0 . 

i-l  s.  1 1 
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Now  u € V and  thus,  by  Theorem  2.3,  u is  tho  limit  in  W of  functions 

oo  «}i  2 h 

Vj  « CqCTq).  Consequently,  4j  v ■ 4^  (u  A^u)  is  the  limit  in  W of  tho  functions 
A^Na'A*^)  < co^Sh  proceeding  to  the  limit  we  find  that 


1 / p(D.u)  (D.  A,  v)dx  » 0 , 

i-1  S.  1 1 1 

'o 


so  that  the  first  term  on  the  right  of  (*)  is  roro. 

Thus,  using  the  Cauchy-Schwar*  inequality  and  Lemma  5.4, 


N2  < 2c,  f / p|nAD.u|  | Au | dx  , 


< 2c1(21/’2N)||D1u  i L2(m||  , 


and  hence 


N < 4Clj|u||  . 


This  bound  holds  for  all  sufficiently  small  h.  In  consequence,  appealing  to  Lemma  5.5  and 
(i)  above,  wo  see  that  u,  |s  and  u,  . |s  » u,,_|s  belonu  to  I,‘(S). 

1 i *1  1 « P 


Proof  of  (iii)  i Prom  Tlieorem  S. 8 of  Gilbarg  and  Trudingor  we  know  that  u,  exists  in  S 
and  satisfies 


u’22  “ " u,U  + Ju<2/x2  . 

-1  2 

n.o.  in  S.  If  o « x0  then  u,„  and  u.^/x,,  belong  to  b“(S),  and  thus  so  does  u,,,. 

Proof  of  ( iv) ; It  follows  from  (ii)  and  (iii)  that  u belongs  to  tho  Sobolev  space 
Prom  the  Sobolev  embedding  theorems  u e CIS)  (Adams  11875,  Theorem  5.4)). 

Furthermore,  since  u c V,  it  follows  from  Theorem  2.2  that 


u(xltx0)  > 0 as  x.,  0 
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t 


„ ! 


• i 


j 

;.  i 


I 


for  almost  all  x^.  Since  u is  continuous  on  S we  conclude  that  u ■ 0 on  9S  n rQ. 


proof  of  (v) : This  is  an  immediate  consequence  of  Theorem  2 of  Huber  (1954) . 


O 


Theorem  5.6  informs  us  that  u is  well-behaved  near  r^.  Away  from  r the  operator 
A is  well-behaved.  There  are  many  results  on  the  regularity  of  solutions  of  variational  in- 
equalities for  coercive  operators  (I.ewy  and  Stampacchia  (1969),  Frehso  (1972),  Gerhardt  (1973), 
Uresis  and  Kindorlehror  (1974)).  However,  there  is  a difficulty  to  bo  overcome  before  these 
results  can  bo  applied:  the  function  $ is  not  smooth.  This  is  because  when  wo  integrate 
along  the  characteristics  of  tj>  as  in  Theorem  4.2  wo  find  that  certain  points  in  SI  lie  on 
two  characteristics.  This  is  best  soon  by  considering  It  follows  from  an  analysis  of 

the  characteristic  equations  (4.11)  to  (4.15)  that  if  dR/dx^(Xj)  > 0 then  the  characteristic 
starting  at  (x^,  R(x1>)  intersects  I'22  (Cryor  (1979a)).  On  the  other  hand,  the  character- 
istic starting  at  (l,.  R(L)>  coincides  with  I'22.  At  points  which  lie  on  two  characteristics, 
$ must  be  taken  to  be  the  larger  of  the  two  valuos  obtained  by  integrating  along  the  charac- 
teristics. 

The  motivation  for  the  following  arguments  is  as  follows.  We  cannot  prove  directly  that 
u « R-'^fn)  because  \J>  / ir,t5(0).  However,  the  discontinuities  of  occur  in  the  upper 
right  corner  of  fl  where  in  gonoral  the  material  is  elastic  and  u > <!>.  We  therefore  seek 
to  replace  $ by  a smooth  function  ^ which  agrees  with  t!1  when  u • $. 

hot  u^  denote  the  solution  of  the  elastic  problem  corresponding  to  the  elastic-plastic 

problem.  That  is,  Au^  • 0 in  0 and  u^  satisfies  the  boundary  conditions  (1.7)  to  (1.9). 

u satisfies  u < K and 
e e 


whore 


a(u  , v -u  ) » 0 for  v c K 
eeo  e o 


K “ (v  t V : v » T/2n  on  in  the  sense  of  H^tsl)) 

o o e 1 


(5.12) 


(5.13) 


As  in  Theorem  5.6  wo  conclude  that  u is  smooth  in  some  S.  . From  the  standard  thooiv  of 

e h 

elliptic  equations  we  conclude  that  is  at  least  twice  continuously  differentiable  vn 

fl.  , so  that  u . c“  (i)) . 
h e 
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Lemma  5.7. 


r- - 


u > u 
— e 


ft- 


Proof:  See  Stampacchia  11965)  and  Cryer  and  Dempster  [1978].  Let  C «■  max(u,u  ).  The 

© 

theorem  will  be  true  if  we  can  prove  that  u » c. 

Now  (as  in  the  proof  of  Theorem  j.2)  ; ( K and  since  u satisfies  (4.7), 

a(u-c,  u-C) 

- a(u,  u-C)  + a(C.  C-u)  , 

< + a(C , C-u)  , 

■ a(u  . C-u)  + a(u  -c,  u-C)  . 
e e 


But, 


where 


a(u  , C-u)  * a(u  , v -u  ) * 0 , 

e e e o 


V * u + C-u  £ K 
e e e 


Furthermore,  either  C ■ u or  C ■ ue  and  so 


a(uo-C,  u-C)  “ 0 . 


Thus,  a (u-C,  u-C)  < 0 and  we  conclude  that  u ■ C. 


L) 


Definition  5. 1. 


2,® 


ij,  satisfies  Condition  C if  there  exists  ij>  £ H ' (0)  such  that  -A  « i>>  whenever 


u < ib  and  ib  < ili  < u whenever  u > ft. 
G — v y — *c  — e e r 


Remark  5.7. 


Condition  C can  be  checked  knowing  only  ij>  and  ue,  both  of  which  can  be  evaluated 
fairly  easily  and  do  not  depend  upon  u.  Conditibn  C is  satisfied  in  some  practical  cases 
(Cryer  [1979a]). 

If  ijt  satisfies  Condition  C we  introduce  the  variational  inequality  with  unique  solution 

u : Find  u £ K such  that 
c c c 


a(u  , v -u  ) > 0,  for  v c K , 
c c c — c c 


(5.14) 
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if  v *■  '■  * • *9%  p-’WW  J 


K = {v  e K : v > Ip  } 
c c e c — c 


(5.15) 


Lemma  5.8. 

If  t(i  satisfies  Condition  C then  u = uc> 

Proofs  u and  u are  the  unique  solutions  of  the  variational  inequalities  (4.7)  and  (5.x4), 
respectively. 


Let  v e K . Then  v e K since  v > <6  > ip.  Thus, 

c c c c — C “■ 


a(u,  vc~u)  >_  0 if  vc  £ Kc 


Furthermore,  u >_  ip  because  ei ther  ug  <_  ip  in  which  case  u >_  t(>  = <J<c  or  ue  > <p  in 
which  case  u > u * i)/  . 

0 C 

We  conclude  that  u also  solves  the  variational  inequality  (5.14) , so  that  u = u . 


Theorem  5.9. 

If  ip  satisfies  Condition  C then  u ■ u e H2,p(fl)  n c (S.  ,,) , for  any  e (1,»). 

c n0/<: 

I 00 

Outline  of  Proof;  It  was  shown  in  Theorem  5.6  that  u|S  e C (S)  where 
S = {x  = (x^.Xj)  s L/4  < x^  < 3L/4,  0 < x2  < h0/2}  . 

The  restrictions  on  the  length  of  S can  be  easily  removed  by  enlarging  0 to  0 by  re- 

00 

flection  as  in  Lemma  5.1.  We  conclude  that  u e C (S,  ,_) . By  Lemma  5.8,  u = u . 

V2 

Now  let 

K,  = {v,  c H*(fi.  ..)  : v = u = u in  the  sense  of  H1  ( Jl)  on 

1 1 "o'2  ^ C 

3SU  n 30  ; v.  = T/2rr  on  T,  in  the  sense  of  H1(S)),-v  > i p } , 

h q/2  hg/2  1 1 i — c 

and  let  a,  be  defined  on  H1  (il.  ..)  x H*(fl,  ,„)  by 
1 h0/2  h0/2 


al(vl'Wl!  = / p Div1Diw1dx 

V2 


The  variational  inequality:  Find  u^  satis: ying 


al^ul'  vl"ul^  - °'  for  a11  V1  K1  ' 

has  a unique  solution  u^. 

Now  let  u = u|fl  -2.  Then,  u^  £ K^.  Furthermore,  for  v £ Kj^  let  v £ K be  obtained 
o'*  * L 

by  setting  v(x)  = u(x)  for  x £ S.  ...  Then 

no'* 

a1(u,  v-u)  « a(u,  v-i.)  £ 0 , 

so  that  u also  solves  the  variational  inequality  (*) . Thus  u * u|ft  ■ u.  . 

Since  uls  ..  is  smooth,  it  only  remains  to  show  that  the  solution  u of  (*)  belongs 

V2  1 

to  H2,p(0).  "rom  Condition  C we  know  that  i(>c  e H1  (0)  n H ,p(fl).  Furthermore,  a^  is  well- 

behaved  on  ft  ...  The  regularity  results  in  the  literature  (Lewy  and  Stampacchia  [1969, 
ho/2 

Theorem  3.1),  Brezis  and  Stampacchia  (1968,  Corollary  II. 3),  Stampacchia  [1973,  Theorem  6.4)) 

are  not  immediately  applicable  because  they  consider  the  case  . It  is,  however, 

2 . 9 

clear  that  the  arguments  can  be  modified  so  as  to  show  that  u,  £ H ,p(fi.  . 

1 h /2 


,JS*  V 


6.  Bounds  for  grad  u. 


We  introduce  tho  elliptic  operator 


Mv  « div(x2  grad  v) 


Lemma  6.1. 


Lot  v c n satisfy  Av  ■ 0 in  Let  w * v-l. 

X2 


sup  w « sup  w 
«h  3nh 


Proof i Since  av  *«  0 in  (l  wo  havo,  using  summation  notation, 


3 v'ii  4 v,2 

X2  X2 


so  that,  by  differentiation, 


V'iij  " <^V'2)'j 


I 2 [v,,  V,,] 


Then  we  C2(n.  ) n c°(n.  ) . Also, 
n n 


Mw  » (x2(v,^vf:j/x2)  , (i  , 

3 4 

- (2v, .v, , ,/x  ) ,.  - 4(v,  v,,/x  ),  , 

j ij  2 l 3 3 2 2 

- 2<v,i.r/x3  V av.jV.uj/x’  - 


- 6v,.v,2j/x42  - 0v(jv,2j/x4  ♦ 


+ 16v' jV, j/x2 


Using  (6.3)  to  roplaco  v ...  and  collecting  terms, 

> tr  J 
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Mw  - 2 "YXj  + 6v'jv'2j^x2  “ 

- 6(v,2)2/)<2  - 6v>jv»2j'/x2  “ 8v'jv'2j'/X2  + 

+ 16{v^)2/Kj  . 

“ 2(v,i;j)2/x2  - - 

- 6(v,2)2/x2  + 16(vtj)2/x^  . 

Thus, 

Mw  “ 2(v'lj  > 2/X2  + 2fv'2)2/>x2  + 8(v'x)2//x2  + 

2 

+ 2x2lV'2j  /x2  " 2v'j/X2J  ' 

> 0 . (6.5) 

Sinco  M is  an  olllptic  operator  in  S)^  wo  conclude  from  tho  maximum  principle  (Cilbarg 

and  Trudingor  [1977,  p.  31])  that  w attains  its  maximum  on  tho  boundtwy  of  S)^. 

□ 

Thoorum  6.2. 

If  k s k^<  if  t natiuflos  Condition  C,  and  if  u is  tho  solution  of  tho  variational 
inequality  1-1.7)  thon 

w - l^rft2  < k,  in  SI  . (6.6) 

x., 

o p 1 — 

Proof;  Wo  know  from  Theorem  5.9  that  u t H*'1  (SI).  In  particular  u i C (11).  l.ot 

S)^  » {x  c ft  : u > *M  , 

H0  « {x  <s  SI  : u ■ (Ji) 

In  S1Q  wo  have  that  grad  u « grad  i)>  so  that,  from  tho  definition  of  \f> 

w « k,  in  Sip  . (€>./) 

4 

Prom  Theorem  5.6  wo  know  that  u «•  x0v  in  . Thus,  for  some  hj  h^  wo  have 
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w < k/2,  in  S 

hl 

Now  consider  the  set  (!  = Q ,/S,  . Applying  Lemma  6.1  we  conclude  that 

1 

max  w = max  w . 

n an 

Now,  3n  consists  of  several  components  which  we  consider  in  turn. 

(i)  P.  “ an  A as,  . Using  (6.8)  we  conclude  that  w < k on  ?„. 

o hj  0 

(ii)  r2  = an  ft  r2.  Since  we  can  use  Theorem  5.1  to  enlarge  n,  we  know  that 

max  w < max  w . 


(iii)  r+  «■  an  n an  . Since  w is  continuous,  it  follows  from  (6.7)  that 

max  w = k . 
f+ 

(iv)  = an  n Pj^.  On  r wo  have  u ■ <Ji  ■ T/2 it.  Since  Au  ■ 0 in  fi,  and  u <_  T/2ti 

on  an,  it  follows  from  the  maximum  principle  applied  to  Au  = 0 that  u T/2ti  on 
n.  Now  consider  a point  on  f^. 

Let  t and  n denote  the  normal  and  tangential  directions,  so  that 


| grad  u|2  = u2  + u2  . 

n t 


Since  u *>  T/2w  on  P^,  wo  have  u^  = 0.  On  the  other  hand,  along  the  inward  normal  n we 


<_  u £ T/2n  . 


Remembering  that  ij>  = T/2n  on  , we  conclude  that 


1 9 tad  u|  = |un|2  < |i|<nl  = kx2  , 


so  that  w £ k on  . 


Remark  6.1. 


Theorem  6.1  is  analogous  to  the  result  of  Brezis  and  Sibony  [1971!  for  the  elastic 
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